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Abstract 

Quantum computing takes fully advantage of the superposition principle to 
increase greatly (even exponentially) the speed of calculations, relative to the clas- 
sical approach. The Deutsch-Jozsa algorithm is the simplest quantum algorithm 
illustrating this power. Unfortunately, the standard derivation involves several 
ingenious steps which usually leave students feeling that they could never have 
figured out the algorithm by themselves. We present here a different formulation 
of the problem which allows students to derive the algorithm using only basic 
knowledge of quantum mechanics and linear algebra. 
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1 Introduction 



Some of the most exciting applications of modern quantum mechanics are in the field 
of quantum computing. As such it would be worthwhile to introduce the topic in un- 
dergraduate quantum mechanics classes, even only briefly. Although the basic idea of 
quantum computing is simple, to harness the power of the superposition principle to 
speed up calculations, the actual design of quantum algorithms is far from intuitive. 

The simplest and among the most well-known is the Deutsch-Jozsa algorithm which 
we will describe in details below (introduced in Refs.jT] and [2] and subsequently im- 
proved in Ref. [3]) . The basic idea of the algorithm is to show how, for a specific type of 
operator, treating the system quantum-mechanically and making a clever choice for the 
input quantum state allows one to gain information in fewer operations than would be 
possible with any classical system. This is of course the essence of the power of quantum 
computing but the beauty of the Deutsch-Jozsa algorithm is that it exhibits this power 
in the simple setting of a two dimensional Hilbert space and without requiring any other 
mathematical concept than basic linear algebra (as opposed to Shor's algorithm, for 
example, which involves quantum Fourier transforms). As such it is the ideal example 
of quantum computing for an introductory quantum mechanics class. 

It is fairly straightforward to verify that the particular quantum superposition pro- 
vided by the algorithm does indeed allow to gain information in fewer steps than would 
be possible classically and checking this, in itself, can be assigned as a homework prob- 
lem. However, it is not obvious how one could have seen that this input quantum state 
is the right choice to start with. In other words, it is easy to check that the solution 
works but it is hard to see how one could have figured it out in the first place. After 
checking that the algorithm works, students therefore often ask: How did anyone ever 
come up with this solution? Of course, the discoverers of the algorithm had a lot of ex- 
perience in quantum mechanics and are very clever, which is why they were able to find 
the solution. Indeed, even the standard presentation of the solution (presented below) 
involves several steps which are nontrivial and which usually leave the students feeling 
that they could never have thought of that by themselves. 

From a pedagogical point of view it would be more interesting to have the students 
figure out the algorithm by themselves as opposed to simply applying it. The goal of 
this short article is to show how the problem can be posed in such a way that it can be 
solved by students taking an introductory quantum mechanics class and using none of 
the clever tricks involved in the standard presentation but only basic linear algebra. 

In the first section we review the problem to be solved. In the following section we 
give the standard presentation of the Deutsch-Jozsa algorithm. In the last section we 
present a way to state the problem such that the algorithm can be "discovered" by the 
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students using only simple algebra and a basic understanding of quantum mechanics. 
The solution is then presented. 



2 The problem to be solved 

Consider for now a classical bit whose states we will denoted by or 1. In this section we 
consider only classical systems. Consider a function which acts on a single bit [I]. Let's 
define a constant function one which always returns the same value no matter what the 
input is. Clearly, there are two such constant functions: the one which always returns 
0, which we will denote Cj, and the one which always returns 1, denoted by Cjj. 
In other words, we have that the action of Cj on a bit is given by 

C/(0) = 0, (2.1) 

Cj(l) = (2.2) 

whereas 

C n (0) = 1, (2.3) 

C H (1) = 1. (2.4) 

Now we consider the so-called balanced functions, which are balanced in the sense 
that they may return or 1 depending on their input. There are also two possible 
balanced functions. We will call the first one Bj, with action 



5/(0) = 1, (2.5) 
= 0. (2.6) 

The second balanced function, Bjj is simply the identity operator: 

Bjj(O) = 0, (2.7) 
B H (1) = 1. (2.8) 

Clearly, the balanced operators are invertible; given one of the two operator and a 
certain output one can determine uniquely what the input was. This is obviously not 
the case for the constant functions. 

Another operation that is required before we can proceed is the exclusive OR (XOR) 
operation © which is simply the addition of two bits (in base 2): 

0©0 = 1©1 = 0©1 = 1©0 = 1. (2.9) 
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Another way to think of XOR is as a parity operation. Indeed, if one assigns a positive 
parity to and a negative parity to 1, the operation a© b essentially computes the parity 
of the product ab. 

Consider now an operator that takes two bits x and y as an input. The first bit, x, 
is used as a control bit and is therefore left unchanged. The second bit is replaced by 

y^f{x)®y (2.10) 

where the function / is one of the four functions described above. We will represent such 
an operator by the symbol Tf where / is either Bi, Bjj, Cj or Cn- We can therefore 
represent the operation of Tf on two bits by 

Tf (x,y) = (x,f(x)®y). (2.11) 



Note that the operators Tf are invertible; they are actually their own inverse. Indeed, 
if we apply any which one twice in a row we get 

T(T(x, y)) = T(x, f{x) ®y) = (x, f(x) © f(x) © y) = (x, y) (2.12) 

where we have used that for any bit x we have f(x) © f(x) — and therefore f(x) © 
f(x) ®y = y. 

To consider a specific example, pick an operator with the function / being the first 
balanced function Bj. If we feed to this operator the two classical bits x — 1, y — 1, the 
output will then be 

T Bl (l, 1) = (1,1© 5,(1)) = (1, 1 © 0) = (1, 1) (2.13) 

Now the problem to be solved can finally be stated. Imagine being provided one of 
the operators Tf but not being told which of the four functions it uses to calculate the 
output. Such an operator with / being unspecified will be called an oracle and will be 
represented by Of with the understanding that / is unknown. You may only feed it 
pairs of bits and examine the result coming out. The question to answer is the following: 
what is the minimum number of times the oracle must be run in order to determine with 
certainty if the function it uses is balanced or constant? 

Of course, another obvious question would be to ask how many trials would be 
required in order to tell which of the four functions is being used, and we will get 
back to this point below, but this is not the question addressed by the Deutsch-Jozsa 
algorithm. 
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It is not hard to convince oneself that two trials are required in order to being able 
to tell if the function is constant or balanced. This will be obvious if we get back to our 
above explicit example where the input was (1,1) and the output was also (1,1). In our 
example we computed the output using f = Bj but let's pretend that we are only given 
the input and the output and that we don't know what the function / is. What can 
we tell about / from this single run? The first bit is always left unchanged so we can't 
learn anything from it. The only information we have is from the second bit. We know 
that the function used by the oracle satisfies 

/(l) + 1 = 1 (2.14) 

which implies 

/(1) = 0. (2.15) 

This is the only information provided by the above test run. But there are two functions 
having this property: Bj and Cj. Therefore, given only that the input was (1,1) and the 
output (1,1), we can only say that the function is either Bj or Cj. To determine which 
is actually used, we would have to run a second trial. For example, feeding the state 
(0,1) will produce the output [5] (0,0) if the function is Bj as in our above example but 
would have produced the output (0, 1) if the function had been Cj. 

It is clear that running the oracle only once for some input (x, y) gives either the 
information f(x) = 1 or f(x) = 0. In either case, there are always two possible functions 
to choose from, one balanced and one constant. The oracle must be run a second time 
to pick which one. 

The above conclusion is inescapable if the bits are classical. The magic of quantum 
computing is that if we allow the oracle to act on quantum bits (linear superposition of 
classical bits) it is possible to determine if the function is balanced or constant with a 
single run of the oracle (but one can't determine which of the two balanced or constant 
functions is being used). This can be done if one is clever about the choice of the two 
qubits being fed to the oracle. The determination of this clever choice of the input 
quantum state is what the Deutsch-Jozsa algorithm accomplishes. In the following 
section we review the standard presentation of the algorithm. 

3 The Deutsch-Jozsa algorithm: standard presenta- 
tion 

Te Deutsch-Jozsa algorithm provides a clever choice of the two qubits to input in the 
oracle in order to determine if the function is constant or balanced with a single run. 
Instead of simply stating the answer we will try to provide the motivation for this choice. 
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This time the oracle takes as input a tensor product of two qubits and outputs 
another tensor product of two qubits. The effect of the oracle on a tensor product of 
the basis states is the following: 

O/(|0)®|0» = |0)®|/(0)©0>, 

0/(10)0 11)) = |0> ® 1/(0) © 1), 

O/(|l>®|0» = |1)®|/(1)©0), 

0/(11)0 11)) = |1)®|/(1)©1). (3.16) 
A general input state may be written in the form \a) © |6) with 

| a) = ai|0) + a 2 |l) 

|6> = 6x|0) + & 2 |l). (3.17) 

Since the oracle is a linear operator, we immediately obtain the result for the appli- 
cation to an arbitrary state to be (using Eq. (13.171) for the two kets) 

O f (\a)®\b}) = a x b x 0/(|O)® |0» + ai & 2 O/(|0) ® |1» 

+ a 2 b x O f (\l) © |0)) + a 2 b 2 O f (\l) ® |1)) 

= ai6i|0) © |/(0) © 0) + ai& 2 |0) © |/(0) © 1) 

+ a 2 h\l)® |/(1)© 0) + a 2 6 2 |l)® |/(1)© 1). (3.18) 

Note that in general this is not a tensor product of the control qubit \a) times some 
other ket, i.e. it cannot be written in the form \a) © |c) for some ket |c). We therefore 
see that in general the control qubit \a) gets entangled with the second input qubit. 

The first trick is to observe that one way to distinguish balanced and constant func- 
tions is to compute /(0) © /(l). Indeed, the result is 

/(0)©/(l)=0 (3.19) 

if the function is constant and 

/(0)©/(l) = l (3.20) 

if the function is balanced, as can easily be verified. 

At first, it seems as if using this will require running the oracle twice since the 
function apparently has to be applied twice (once on and once on 1). But this is not 
necessary if one works with quantum bits at the condition that the function appears 
through a phase, as we will now show. 
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Schematically, we are trying to have a situation where a relative phase of the desired 
form will be generated between the two input bits |0) and |1). In other words we would 
like to find a qubit \b) which has the properties 

O/(|0)®|6)) = (-1) /(0) |0>® 1 6) (3.21) 

and 

O f {\l)®\b)) = (-lp)|l>® \b). (3.22) 

As we will see below, it is not difficult to find a ket \b) that has those properties. 
But before presenting the (almost obvious) answer, let us first discuss why having such 
a state allows one to distinguish constant from balanced functions in only one step. 

Assuming that we have a state \b) satisfying the above properties, applying the oracle 
to the product \a) ® \b) where \a) is an arbitrary qubit will give 

O f (\a)®\b)) = oi O f (\0)®\b))+ a 2 O f (\l)®\b)) (3.23) 
= (-l) /(0) ai \0)®\b) + (-l) m 02 1 1)® 1 6) (3.24) 



(_l)/(o) UlO) <g> \b) + (-l)-m+fV) a 2 \l) ® \b)j (3.25) 

( ai \o) + (-i)m+m a2 \i)) \b) (3.26) 



where have used (— 1)~^ ) = (— As usual, the overall phase is unimportant. 

Therefore, using that /(O) © f(l) = 1 if the function is balanced (see Eq. fl3.20p ). the 
final result is simply (after discarding the overall phase) 



O f (\a)®\b)) = f ai |0> - aall) 1 ® \b) 

= ai|0) ® 1 6) — a 2 \l) ® 1 6) (for a balanced function). (3.27) 

and if the function is constant, using /(0) + /(l) = 0, we simply get 

O f (\a)®\b)) = \a)®\b) 

= ai|0) ® 1 6) + a 2 1 1) ® 1 6) (for a constant function.) (3.28) 

If we impose that the two possible outcomes are orthogonal so that they can be 
distinguished, we find that we must have |ai| 2 — \a 2 \ 2 = which implies that ai and a 2 
may only differ by a phase: 

a 2 = e*% (3.29) 
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Therefore, the normalized control qubit must be 



|0) + e* e |l) 

| a) = -j= (3.30) 

and the input state that must be fed to the oracle is 

|n\ + e ie \l) 

\a) ® |6> = |U; ^ 1 ; ® \b) (3.31) 

with the ket \b) having the properties (I3.2ip and (I3.22p . 

We have accomplished the goal we had set for ourselves: we have found an input 
state such that running the oracle only once will determine if the function is constant 
or balanced. Indeed, if the we use the above state as input, we simply have to project 
the output on the bra (a| ® (6|. If the result is one, the function used by the oracle is 
constant. If the result is zero, the function is balanced. 

All this relies of course on finding a state \b) satisfying the properties (I3.2ip and 
(I3.22p . It is not difficult to guess what the answer is: 

I*) = (3.32) 
Consider the case /(0) = 1 (so the function is either Bj or Cu ). Then we have 

|0)<8> |/(0)®0) |0) <8) |/(0) ® 1) 
|0)_g_|l) |0)_gj0) 

= -|0)®|6) (3.33) 

which we may write as (-1) /(0) |0) <g> |6). 

It is easy to verify that if /(0) = (so the function is either Cj or Bjj), the state 
|0) ® \b) is left unchanged by applying Of so that we may still write the result as 
(— |0) <8> \b). The conclusion is therefore that for any of our four operators, we 
obtain 

O f (\0)®\b)) = (-l) m \0)®\b) (3.34) 
with | b) chosen as in Eq. (jX 
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It is easy to check also that for any operator, one finds 



/ (|l)®|6)) = (-l)/( 1 )|l)(g)|6). (3.35) 

So the state (13.321) satisfies the relations (I3.2ip and (I3.22p . 

Our final result is therefore that the input we must feed to the oracle is 

|0> + e »|l>W|0>-H>\ (3.86, 



V% J V y/2 

With this choice of input, the outputs corresponding to an oracle using a balanced 
function will be orthogonal to the output corresponding to an oracle using a constant 
function. As explained above, one can determine if the function is constant or balanced 
by projecting on the bra (a\ <8> (b\. 



4 An elementary derivation 



The presentation given in the previous section emphasizes the key ideas on which the 
algorithm is based but students (and non-students alike!) may feel that it relies on 
too many clever guesses. We will present now the problem in such a way that it can 
be solved "mechanically" using only basic linear algebra and none of the clever tricks 
needed in the previous derivation. It is of course more algebra intensive and may be 
considered less satisfying from a conceptual point of view by some. But it has the 
pedagogical advantage that it can be solved by students with only a basic knowledge 
of quantum mechanics who will therefore have the satisfaction of having "rediscovered" 
the algorithm by themselves. 

We will represent the tensor products of two qubits as four component column vectors 
using the convention 



|0> ® |0) 



|1) <8> |0) 



/ 1 \ 





V o y 
/o\ 



1 



|0) <g> |1) 



II) ® II) 



1 



v o y 

/o\ 




V 1 / 



(4.37) 



(4.38) 
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The first step is to represent the four possible operators Ff as four by four matrices. 
It is now easy to construct explicit representations of those operators since we know how 
they act on each of the four basis states. For example, 



F ClI (\0) ® |0» = |0> ® \C n (0) © 0) = |0> ® |1) 



(4.39) 



and so on. One finds 



F Cl 



F Br 



( 1 ^ 

1 

10 

\ 1 J 

( o 
1 

1 



F CI 



/o 
1 



1 \ 




1 \ 








\ 1 J 



1 
\0 1 0/ 

/ 1 \ 



F Bj 



1 





\ o o i 





1 



(4.40) 



(4.41) 



0/ 



Now we are ready to state the problem in such a way that students in an introductory 
quantum mechanics class could solve without any subtle trick, only simple algebra, and 
therefore recover the solution provided by the Deutsch-Jozsa algorithm. 

Let's consider an arbitrary input state 



( ci N * 

c 2 

\ C4 y 



(4.42) 



with 



M 2 + 



c 2 | 2 + 



c 3 | 2 + 



|c 4 f = I. (4.43) 

The goal is to choose the coefficients in such a way that, given an oracle, we will 
be able to tell with only one application whether the function used by the operator is 
balanced or constant. The way to obtain this is clear: we must choose the coefficients 
of our input state such that if we run it through each of the four operators, the output 
of each balanced operator must be orthogonal to the output of both constant operators. 
This gives four conditions which we can schematically write as 



{FcMObA) = 0, 
(FcMOb u ^) = 0, 
{F Ci MOb^) = o, 



(4.44) 
(4.45) 
(4.46) 
(4.47) 
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This is all there is to it. If one can find coefficients satisfying these four conditions (and 
of course we know that there is at least one solution), the problem has been solved. 
Equation (I4.44p yields the following condition: 

C 2 C 1 c l c 2 "I" | c 3| 2 "I" | c 4| 2 

= 2Re (dc* 2 ) + |c 3 | 2 + |c 4 | 2 = (4.48) 

whereas equation (I4.45P gives 

2Re(c^c 4 ) + |ci| 2 + |c 2 | 2 = 0. (4.49) 

The third and fourth equations, (14.461) and (14.47}) . end up being exactly the same as the 
two above. So we need to solve those two equations with the constraint that the state 
be normalized, equation (l4.43p . To simplify the problem, one may at first try to see if 
there exists a solution with all coefficients real. Our set of equations is therefore 

2 Cl c 2 + c 2 + c 2 = 0, (4.50) 
2c 3 c 4 + c 2 + c 2 = 0, (4.51) 
cl + c 2 2 + 4 + cl = 1. (4.52) 

If we isolate c| + c\ from equation (l4.50l) and use this in equation (I4.52[) . we obtain 

c\ + c 2 - 2c lC2 = 1 (d - c 2 ) 2 = 1 (4.53) 

which gives 

ci = c 2 ± I. (4.54) 
Similarly, if we isolate c\ + c 2 from (14.511) and put this in equation (14.52j) . we get 

c 3 = c 4 ±l. (4.55) 

4.1 First case: c\ = C2 + 1 and C3 = C4 + 1 

Let's first pick the cases 

ci = c 2 + l, (4.56) 

c 3 = C4 + I. (4.57) 

If we plug those values in the normalization condition ( 14.521) , we get a quadratic formula 
for c 2 in terms of C4: 
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If we impose that c 2 be real, the discriminant must be larger of equal to zero: 

- 4c 4 6 2 - 4c 4 - 1 > (4.59) 

but it's trivial to show that the only solution is for a strict equality (in other words, this 
describes a parabola with its apex on the horizontal axis) which is solved by 

c 4 = -~ (4.60) 

Putting this back in equations (I4.55p . (14.561) . (I4.57P and (14.581) . we final get 

ci = -c 2 = c 3 = -c 4 = - (4.61) 

as our first solution. Written as a tensor product of two kets, this corresponds to 

|(|0>«> 10) ~ 10) ® 11) + 11) ® 10) - II) ® |1)) = ( |0) ^ 1)} ) ® ( |0) ^ 1)} ) (4-62) 
which is indeed a special case of equation (l3.36l) with 6 = 0. 

4.2 Second case: c\ = C2 + 1 and C3 = C4 — 1 

In that case, one finds obviously the same quadratic equation as before except that with 
the opposite sign for the term linear in c 4 : 



^ = -l±V-4cI + 4c 4 -l (4 63) 

Again, imposing c 2 to be real leads to solving 

_ Ac\ + 4c 4 - 1 > (4.64) 

which, again, turns out to have a unique solution, corresponding to a strict equality. 
The result this time is c 4 = \ which yields 

ci = -c 2 = -c 3 = c 4 = - (4.65) 

giving the state 

\ U) ® |0) - |0) 8. [1) - jl) ® 10) + [1) ® = ( |0) ® ( |0) ^ 1)} ) (4.66) 

which is indeed the other solution of equation fl3.36l) with real coefficients, i.e. 9 = n. 
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4.3 The other two cases 



We still have to consider the cases c\ = c 2 — 1 and C3 = C4 ± 1. As one would expect, 
these two cases give the same two results we obtained above excepted for an overall 
irrelevant factor of — 1 . 

5 A Bonus 

One obvious question that comes to mind when learning about the Deutsch-Jozsa algo- 
rithm is whether it would be possible to not only determine if the oracle uses a constant 
or balanced function in only one step but in addition to determine uniquely which of the 
four functions is used. Since the standard presentation relies on several clever guesses, 
one might wonder if more clever tricks could be used to identify uniquely the function. 
The answer is easy to find using our "mechanical" approach. The question becomes 
whether it is possible in addition to the conditions (14.441) - (I4.47P to fulfill the extra 
constraints 

(TcMFc u ^) = (5.67) 
(FbMFb u ^) = (5.68) 

which both give (going back to complex coefficients for now in order to be more general): 

c*c 2 + c* 2 c\ + c%cn + C4C3 = (5.69) 

which may be written as 

Re(c*c 2 ) + Re(c*c 4 ) = 0. (5.70) 

But this condition is clearly inconsistent with equations (I4.48P and (14.491) which com- 
bined imply 

Re(c*c 2 ) + Re(c*c 4 ) = --. (5.71) 

Therefore we see that it would be impossible to determine uniquely the function used 
by the oracle in a single step. 

6 Conclusion 

We have shown how it is possible to formulate the problem solved by the Deutsch-Jozsa 
algorithm in such a way that the solution can be found using only basic concepts of 
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quantum mechanics and simple linear algebra. This brings the rediscovery of the algo- 
rithm within the reach of students taking an introductory class in quantum mechanics 
as part of an assignment set, for example. 
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